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ABSTRACT 


Using a computer program a partly optimized 
configuration for a superconducting version of side 
and lift coil system of NASA-MIT prototype is 
presented. Cable size for the mentioned coils 
and also for superconducting drag and magnetizing 
coils regarding the overall computed field is 
determined. 
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Chapter . I 


INTRODUCTION 
1.1 Magnetic Balances 

Magnetic suspension of models in wind tunnels 
is more advantageous than conventional support of 
models for two reasons: first, it has no effect 

on the flow and secondly, it provides fhe capability 
to perform dynamic stability tests. 

Since magnetic suspension of models in large 

wind tunnels with high dynamic pressure, such as 
2 

NASA 1 ft pilot cryogenic # wind tunnel, requires 
high fields, using a superconducting magnetic 
balance would reduce the required power. 

t 

Investigation of using superconducting coils 
for magnetic suspension and balance systems has 
been done at the M.I.T. Aerophysics Laboratory in 
three parts: first, Kraemer (1) studied low A.C. 

losses and found the type of superconducting wire 
that had acceptably low losses. Secondly, Prey (2) 
performed some experiments on A.C. losses in inter- 
acting superconducting magnetic coils. This thesis 
tries to design a superconducting coil array for 
present side force, lift, drag and magnetizing coils 
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of the NASA prototype magnetic suspension and balance 
system. So, it provides a preview of full design 
procedure for such a system. 

1.2 Superconductivity 

It is well known that superconductivity is the 
absence of electrical resistance due to interaction 
between electron pairs and the crystal lattice of 
superconducting material. 

Superconductivity can be extinguished by 
disturbing the electron pairs, such as increasing 
temperature, which results in lattice vibration or 
by an external magnetic field which causes internal 
voltage difference and accelerates the electrons. 

There are two types of superconductors: 

Type I, pure materials, Which have zero resistance to 
direct current but cannot retain superconductivity 
in the presence of an external magnetic field. 

Type II, superconductors, metallic compounds and alloys, 
are more suitable for M.S.B. systems because they can 
tolerate magnetic fields to some extent. 

1.3 The Problem at Hand - 

It was shown in (1) that a suitable conductor for 
an M.S.B. system capable of dynamic suspension of 
models is a cable of fully transposed fine copper wires. 
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each having a single superconducting core and all 
wires being electrically insulated from each other. 

The interacting coil losses using the mentioned 
kind of superconductor was studied in (2) and it was 
shown that the interacting losses are about 10-15% of 
self field coil losses. 

The basic approach to the design of a super- 
conducting coil array for the M.S.B. that has been 
followed in this work consists of these steps: 

First, developing a computer program to compute the 
field due to different coil geometry (straight line 
or circular current element) and then using this to 
find the best coil configuration, which produces the 
desired field. This minimizes the superconducting 
volume and as a result, reduces the overall losses. 
Secondly, computing the field at several wire loca- 
tions inside the coils to determine the minimum 
superconductor cable size. 
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Chapter II 
THEORY 

2.1 Basic Properties 

A basic property of superconductivity which has 
direct benefit in engineering is the absence of 
electrical resistance to direct current. A more 
important aspect is the Meissener effect, which is 
the main difference between a superconductor and a 
perfect conductor. Figure 1 shows this effect; 
after increasing the field above H c , the field in- 
side the superconductor jumps to H c with the loss of 
superconductivity. Upon decreasing the external 
field the inside field again returns to zero upon 
reaching H c no matter what the history of the 
process was. In this regard superconductivity is a 
reversible phenomenon. In contrast to this is a 
perfect conductor, which retains its inside field 
after passing H c and decreasing the field; i.e., a 
superconductor expells all fields as long as it is 
superconducting whereas a perfect conductor could have 
fields frozen in by the onset of perfect condition. 
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Figure 1. Variation of the Inside vs Outside 
Field for a Superconductor and a 
Perfect Conductor, which both have 
a critical field H . 



2 . 2 Superconductors 


a) Type I superconductors have a magnetization 
vs. external field, as shown in Figure 2. There is 
one critical field for them at which both Meissener 
and resistive transition occur. However, super- 
conductivity in Type I is extinguished by a field of 
the order of .1 Tesla or less. So they cannot be 
used in the proposed M.S.B. Lead and Niobium are among 
Type I materials. 

b) Type II superconductors are known with two 
critical fields H ^ and H ^ (Figure 3) . At fields 
less than H ^ there is no flux inside the super- 
conductor. Between H c ^ and H ^ flux penetrates in 
a small depth. The rest of the conductor remains 


in superconducting state. 

If Type II superconductors had no imperfections 
and grain boundaries as low energy sites for this 
penetrated flux (fluxio) , fluxios would move due to 
Lorentz forces and dissipate energy. Fortunately, 
by heat treatment different grain boundaries are pro- 
duced and as a result high pinning strengths are 
attainable. Most interesting engineering super- 
conductors are of Type II. They include metallic 
compounds and alloys, such as NbTi Niobium-Titanium 













alloy, NbZr Niobium-Zirconitjn alloy and Nb-Sn Niobium- 
Tin compound. 

Figure 4 shows the type of relation between 
maximum current density (J c ) and applied magnetic 
field that a Type II superconductor can tolerate 
without transition to the normal state. 

The J c vs. B relation is approximated in two 
regions: a) for small perturbations, in the _mid- 

range of the curve, the simple London approximation 
that J c is independent of B is valid; -b) for lower 
fields the Kim relation is used: 


J 

c 


J 0 + B o 

B + B_ 


t 


J Q and Bq are constants which depend on the flux 
pinning strength of material and self-field of 
conductor. 

2.3 Penetration Depth and Hysteresis Losses 

For a changing field perpendicular to a cylindrical 
wire penetration depth (2) is 



Figure 4. Typical Current Density Relations 
vs. Transverse Magnetic Field for 
a Type II Superconductor Alloy (1) 
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The necessary field change to fu§ly penetrate the 
conductor in the London approximation is 
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Vc a 
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( 1 ) 


For a varying field as (1) shows, since the current 
density is limited, the field that can be shielded is 
also limited. At fields higher than B^, the field 
has reached to the center of conductor and any change 
in external field will appear instantaneously in,- the 
conductor. However, the magnitude of the interior 

field will lag behind an increasing field or a de- 

* 

creasing external field. This causes hysteresis 
losses for Type II superconductor in varying fields. 



Chapter III 


DESIGN 

3.1 The Magnetic Suspension and Balance System 

The present M.S.B. (NASA-MIT Balance) at the 
M.I.T. Aerophysics Laboratory uses three main groups 
of copper-wired coils, cooled with water, with iron 
poles for controlling and producing forces and 
moments in different directions (3) . The three 
groups of coils are as follows: 

a) Helmholtz- coil system (controls B , B ) 

X XX 

b) Saddle-coil system (controls B^ , B z ) 

c) Side and lift force system (controls B , B ) 

. zx yx 

< 

Figure 5 shows the arrangement of these coils. Besides 
these, there is an Electromagnetic Position Sensor Coil 
System, which measures the position of the model for 
controlling it. 

Since the main sources of power consumption are 
coils in groups (a) and (c) and the iron. poles of 
group (c), (Figure 6) , limit the field and makes scaling 
unreliable, it was decided to design coils (a) and (c) 
using superconductors and no iron. 
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The main constraints for such a design are as g 

follows: 

(a) The whole superconducting coils should be 
in one or two warm bore dewars . 

(b) Enough light passage to watch the model 

or use measuring optic devices such as LDA should be 
provided. 

(c) The saddle coil system and position sensor 
coils should remain in their position. 

The remaining space is shown in Figure 7. 

3.2 Computer Program 

The computer program (TABLE) (5,12,13) tabulates the 
magnetic field components due to a series of straight 
line current elements or a series of circular line 
current elements. In the straight line current 
element mode it is limited to the shapes which could 
be approximated by straight lines. Figure 8 shows 
a sample winding that can be handled by TABLE. 

The number of turns in computer program is 
controlled by M and N (number of layers and turns per 
layer in winding cross-section) . The effect of 
changing these parameters and convergence of computer 
results are discussed in Appendix 1. 




Figure 7. Free Space for Design (Unhatched) 
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■rj 3.3 Testing Computer Program 

A model coil similar to the designed lift coil 
system (Figure 9) was made to measure the transverse 
field component on wind axis and compare with TABLE 
predictions for the same coil system (14) . 

Figure 10 shows the results of TABLE together 
with measured results. The computer program over- 
estimates about 10 percent, which could be explained 
by noticing the procedure in which TABLE locates the 
wires in__the winding volume. The cross-section of 
the model coil was not rectangular, so the results 
might be improved to some extent for a rectangular 

t 

winding cross-section. The second important source 
of computer program error is that TABLE approximates 
the curves at winding corners with the intersection 
of two straight lines. 

3.4 Finding the Best Configuration 

The present side and lift coil system with iron 
poles (Figure 6) has a measured performance of: 


Taking I = 385A, the gradient should be 540 gauss/in. 
zx 

This is the basic requirement chosen for each different 
configuration of air core superconducting side and lift 


coils. 



Figure 9. Model Coil for Testing TABLE 


1 Model coil 3 A.C. transformer 

2 Pick-up coil support 4 Pick-up coil voltmeter 















As a first step a irtodel coil with the same 
core form as the iron poles was tested as Model I 
(Figure 11) . Since the performance was not satis- 
factory and winding was difficult, the winding plane 
is not perpendicular to the core form axis so it is 
possible that the wires slip; Model II with the same 
core shape but with the winding plane perpendicular 
to the core axis was tried. Figure 12 shows the 
direction of wires and the different angle of core 
for which computation was done . The disadvantage 
with this model is that due to its acute angle of 
winding with the axis plane, it takes a lot of space 
and leaves the unoccupied space unusable. 

Figure 13 shows Model III. This model came 

* 

about as an approximation to Model II and since it 
uses all available space, it is a better choice in 
most cases. 

Choosing a configuration depends on the 

magnitude of gradient which is needed. For a low 

transverse field gradient (B <200 gauss/in) . 

zx 

Model I and Model II are comparable to Model III but 

at field gradients as high as (B = 500 gauss/in) . 

zx 

Model III is considerably 'superior . For producing - 





X 



Figure 12 . Model II a) Direction of Wires 

b) Different Core Axis 
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500 gauss/in gradient. Model II at 60 core axis 
angle should have 5000 turns in each of 8 coils 
whereas Model III with three-stage optimized coils 
needs 3900 turns in each of 8 locations. This 
should be clear since for extending Model III 
upper stages would be closer to the center than 
upper stages of Model II. 

Model III can be optimized easily if the 
effect of each stage is studied separately. Such 
an investigation shows that at a definite radial 
positioh the highest gradient may be obtained from 
one stage. Table I shows the B /(ampere turn) for 

Z X 

each stage in Figure 13 and a fourth stage that 

might be added. Stage one was discarded after such 

* 

a study because it is far from the optimum radius 
and its contribution to total gradient is negligible. 
Table I shows that the optimum radius is about the 
place where stage 3 is. It seems the optimum radius 
depends on axial distance of two sets besides the 
geometrical parameters of winding. 

It is noticeable that due to symmetry, computa- 
tion of the field of one coil of each stage is 
enough to compute the total field for points on the 
x axis (axis of symmetry). This could save a lot of 


CPU time. 
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Table I 


’ ■ 

Transverse Field Gradient of Different Stages 
of the Side and Lift Coil Model 


Stage 

"Mean 

Radius 

in 

Gradient 

B^auss, 

in 

f 

Ampere Turn 

B /Ampere Turn 
zx 

-• 

1 

■ 

5.5 

15.6 

18x24x260* 

1.39 x 10 -4 

2 

8.5 

102. 

36x24x260* 

4.-54 x 10~ 4 

3 . 

11.5 

190.6 

54x24x260* ' 

5.66 x 10“ 4 

4 

14.5 

189.88 

72x24x260* 

4.23 x 10“ 4 




260A cable current is assumed for design purposes. 


18x24, '36x24, etc. are the each lift coil number of turn; 
i.e. half the total available number of turns. Refer to 
Appendix 2 for more information. 



Comparison of results for the three models shows 
that. Model III is the most advantageous since 

a) It can produce the desired field gradient 

B = 540 gauss/in with the least volume of space 
zx 

and superconductor. 

b) It could be extended for higher field 

gradients. ... . _ 

c) Construction of such a coil system would 
be less troublesome since there is no force in 
windings due to unusual wire direction. 

Figure 14 shows the B component of magnetic 

z 

field for points on x axis. The B gradient at the 

Z X 

* 

center is 586.8 gauss/in with gradient uniformity 12% 
over a sphere of radius 3 in. The average gradient is 
B = 619 gauss/in over a distance from x = -3 to 

ZX 

x = 3 . 

3.5 Magnetizing and Drag Coils 

After dimensions of side and lift coil system 
were found, the magnetizing and drag coil system 
(Helmholtz coils) could be designed according to the 
following considerations: 
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a) Magnetizing and drag coils should produce 
B = 7700 gauss and B = 655 gauss/in according to 

X XX 

the measured performance: 


B 

x 


B 

xx 


20 I 

x 


1.7 


I 

xx 


for present M.S.B. and I = I = 385A 

V V'*' 


b) Drag coils should be inside magnetizing 
for less impedence, so the inner diameter is the 
outer diameter allowed by side and lift coil system. 


c) Magnetizing and drag coils should be 
optimized for the least volume or least inside field, 
since several dimensions are possible trade-off 
between least volume and least field design should be 
effected regarding losses against superconducting 
material expenses. This part should be done after 
loss computation. Note that A.C. losses were not 
computed as a part of this study. 



Computations for these coils were carried ou2 by- 
using a circular current element program (2) . 

Figures 15 and 16 show vs. axial distance for 
magnetizing and drag coils. Magnetizing coils have 
a field component, B =7448 gauss with a uniformity 

X 

of .92 percent over a sphere of radius 3 in and 
corresponding values for the drag coils are B xx = 

634.7 gauss/in at the center with gradient uniformity 

6 . 5 percent over a sphere of radius 3 in . 

Figure 17 shows the overall designed coils with 
present saddle and position sensor coils. All super- 

1 

conducting coils were designed for 260A operation. 

Note that the Helmholtz coils are not separated by 
the theoretically ideal distance. This is also true 
of the NASA prototype Helmholtz coils, which are 
approximately positioned. 

3.6 Magnetic Field inside Coils 

and the Conductor Size 

For finding proper diameter of conductor using 

Figure 4, one should have the magnetic field due to 

all coils inside windings. From these, current density 

should be determined to be less than J . TABLE was 

c 

used for this purpose and the field at different wire 
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Figure 17 legends 


Reference 

Number 


3 

4 

5 

6 


Coil 

Description 

1st stage 
side and lift 

2nd stage 
side and lift 

Magnetizing 

Drag 

Saddle 

EPS 


Req 

No. 


8 

2 

2 


Layer 

36 

52 
38 
5 2~ 


Turn 
per * 
Layer 

42 

42 

50 

75 " ' 


* 

For side and lift force coils these figures are 
half the total of number of turns, accounting 
only for the lift coils t 
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locations was obtained to^ which the wire self field, 
which depends on conductor diameter, should be . 
added later. For these computations maximum 
allowable number of turns were used for neighboring 
coils regarding convergence criteria in Appendix 1. 
Table II includes the maximum field inside each 
stage of side and lift coils, drag and magnetizing 
coils and the percentage contribution of each sfet 
of coils to the overall field. This field is - the 

— ■ 1/2 :r: ‘ 

2 2 2 

magnitude of B; i.e., (B +B +B ) , which is 

x y z 

perpendicular to the wires . The maximum field is 

found from a comprehensive study for each set of 

coils. Figures 18 through 21 show the magnitude of 

the magnetic field inside respective coils. Figure 20 

shows the field distribution in two drag coils., since 

prediction of the worst one is difficult. In. each 

computation absolute value of each field component due 

to side and lift coils were added to the absolute 

value of the same component of the drag and magnetizing 

field. '• • 

*'■ “ + - 

One can define the appropriate conductor size 
either by maximum field in the worst coil, or, by the 



Table II 


Maximum Magnetic Field Inside Coils 
with Percentage of Contribution of each Set 


Coil Description 

% 

Self 

Field' 

S & L Coil 
1st stage 

27.5 

& L Coil 
2nd stage 

35.1 

Drag Coil 

41. 

Magnetizing 

Coil 

34. 

1 


Max 

Cont. of! Cont. of I Mag- J ' Field 

S&LT lstis&L 2nd! netizingj % Drag Tesla 


31.2 


15.4 


11.8 26.8 3.71 


11.1 24.9 


51.1 


3.71 


3.56 


2.81 





































Field magnitude between 
Tesla* 


o 

O 


3.71 - 3.34 
3.34 - 2.97 


Figure 19. 


Field Magnitude in 
Side and Lift Coil 


O 

2nd Stage ' ( 


O 


2.97 - 2.58 

2,58 - 2.21 

2.21 


less than 



Figure 20 





i 


Field Magnitude 
in Drag Coils 


Field magnitude between 
Tesla 


o 

o 

O 

o 


3.71 - 3.34 

3.34 - 2.97 
2.97 - 2.68 
2.68 - 2.21 


less than 2.21 



J 

J 

] 


Field magnitude between 


Sr 

vi 

47 o 

Tesla 

2.81 - 2.55 

a® 

o 

. * • 

2.55 - 2.29 

V 

o 

2.29 - 2.03 


o 

2.03 - 1.77 


C~3 

less than 1.77 
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maximum field in each coil locally. This procedure 
could be repeated for each single coil? i.e., 
dividing each coil into several parts and defining 
the proper conductor size for each segment based on 
maximum field in that segment. The second method 
would be advantageous in reducing superconductor 
volume but would require the purchase of more types 
of cable. 

One can check the initial assumed cable size on 
which the field calculations were based by noticing 
jthat for every cable size trial, the current density in 

superconducting part of the cable should be less than 
the critical current density obtained from a J c vs. B 

relation, such as one in Figure 4. Checking this 
condition leads to a new conductor size and new coil 
dimensions. Therefore, field calculation must be 
repeated until J < J c condition is satisfied. It is 

worth noticing that the field does not change with 
cable size a great deal if one keeps the number of 
turns constant. This may cut down the number of trial 
and error cycles. 



After such an attempt for the proposed configura- 
tion, the proper cable dimensions are presented in 
Table III. The computed maximum fields were increased 
20 percent to take into account the possible computer 
error and other factors. 

The proposed conductor is a cable of solid NbTi 
core filaments, which diameter of filament is .0367 mm 
with copper/superconductor ratio = 1.71/1. This con- 
ductor was tested in (1) and (2) and was proved to be 
the best conductor solution available in market for 
present application. 



Table III 


Final Cable Dimensions 
for Different Coils 


Side and Lift, 1st stage 

Cable _ . 

Dia 

mm 

S SC 2 
mm 

^A/mm^ 

Maxf ield 
x 1.2 
Tesla 

J C 2 

A/mm 

7x6x5 

1.73 

.222 

1170 

4.52 

1250 

Side and Lift, 2nd stage 

7x6x5 

1.73 

.222 

1170 

4.52 

1250 

Drag Coil 

8x5x5 

1.6 9 

• 

.212 

.. 

1229 

4.35 

1338 

. ? 

Magnetizing Coil j 

1 

6x5x5 

1.46 

.159 

1639 

■ 

3.46 

1700 









IV. DISCUSSIONS 


An alternative superconducting version of side 
force, lift, magnetizing and drag coils of the present 
M.S.B. was presented. 

Without any iron poles the side and lift coil 
system requires seventeen times more ampere turns to 
produce the same gradient. However, this ratio could 
be reduced drastically by decreasing the axial dis- 
tance between two parts of side and lift coils. 

The presented design could be more compact if a 
coil form such as one used for saddle coils were used 
for side and lift coils and also if t&e drag and 
magnetizing coils were relocated. 

Choosing conductor « size upon local field (i.e., 
the maximum field in each segment of a coil) might 
reduce the required superconductor volume to some 
extent. By dividing the side and lift coil 1st stage 
into two equal inner and outer parts and choosing the 
appropriate conductor size for each segment (i.e., 
an 8x4x5 cable instead of a 7x6x5 cable for outer part) 
one might save above 16 . 7% in superconducting volume 
for the 1st stage coils. However, it is clear that for 
such a partition all coils of 1st stage should have the 
same field distribution. For example, this procedure 



is not applicable to the drag coils since they have 
a much different field distribution (Figure 20) . 
Perhaps the conductor should be selected for each of 
them separately. 
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Appendix 1 


V 


CONVERGENCE OF COMPUTER RESULTS 


As pointed out earlier, computer program TABLE 
approximates each winding with a series of straight 
line current elements.-. It fills the winding volume 
with the definite number of turns. The choice of 
number of turns depends., on the required accuracy 

and the available computer time. The appropriate 

... re- 
number of turns for every different relative dis- 
tance of coil and the point should be computed 
separately. - 

' t 

Figure 22 shows the results of such a computa- 
tion for a coil having winding parameters (M layers 
and N turns per layer) -' As 2 by 2, 4 by 4 , 8 by 8 
and 16 by 16. (Refer to Figure 8 for the definition 
of M and N.) The vertical axis is relative error to 
the previous step. Figure 22 shows that for such a 
configuration of coil and the point a choice of 16 
by 16 winding results -in an error less than 2 percent. 
This indicates that fdr smaller than 16 wires per layer 
v/ the correct number should be yet but that for larger 
numbers a less costly price value is to substitute an 
equivalent 16 by" 16 array for this particular configura 
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^3 Figure 22. Convergence Characteristic 

, of Computer Program TABLE 
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Appendix 2 

f 

COMPUTER PROGRAM TABLE 
and COMPUTER RESULTS 


Computer program TABLE and computer results for 
the designed lift force, drag and magnetizing coils 
are presented. The lift coils' number of turns is 
half the total available number of turns since the 
second half belongs to side force coils. Therefore, 
the side force field component (B^) distribution is 
similar to lift force component (B ) distribution. 

4i 

In practice, they are powered from a single line 

having the resultant current. 

The complete elliptic integrals in Subroutine 

Circe are computed according to (10). 

< 
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table 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


Program to calculte magnetic field components produced by coils 
consisting of straight line or circular current elements. 

Straight line current elements are counted counter-clockwise about 
the correspondig coordinate di rectio'ns.Al V currents are postive 
counter-clockwise. 


variable name 
im 
jm 
km 
dx 
dy. 
dz 
x{ 1 ) 
yd) 

2 ( 1 ') 
da 
xmu 
co i 1 

xl ,y1 , zl 
x2 ,y2 , z2 
ccur 
i co i 1 
ncoi 1 
ni nt 
mm , nn 
xx, yy, zz 
cl 
xl 
r i 
ro 
1 ay 
nt 1 


Input variable 1 ist 

definition 

number of increments in x direction, 
number of increments in y direction, 
number of increments in z direction. 

"delta x" 

"delta y" 

"delta z" 

x coordinate of starting point for incrementing. 

y coordinate of starting point for incrementing. 

z coordinate of starting point for incrementing. 

demagnetizing constant for the model. 

magnetic permeablity of free space. 

indivdual coil description. ^ 

coordinate of the end points of the straight line 

current elements making up the“doils. 

coi 1 curran t . 

numbe of first coil. 

number of last coil. 

number of interfaces. 

dimenssion of interface array. 

interface endpoints. 

axil distance of solenoid from the center. 

solenoid length. 

solenoid internal dimeter. 

solenoid external dimeter. 

solenoid number of layers. 

solenoid number of turns per layer. 


note: 

In order for coors to work properly, 
ww , wl , w2 ,mus t be dimensioned to at least 

ww( 4 , n i nt ) , w 1 (ntot ) , w2(n tot ) . x x ,yy , zz have the same dimensions 
as ww . 


integer out 

dimension x ( 50 ) , y ( 50 ) , z ( 50) , 

1x1 (873 6) ,y1 ( 8736) ,z1 ( 873 6 ) , x2 ( 8736 ) , y2 ( 8736 ), z2 ( 8736 ) 

2 , conf g (18),curt(500),coi 1(18,28), sum( 12,50),xx(4,10),yy(4,10), 
3zz(4, 1 0) ,ww( 4, 10) ,w1 (8736) ,w2( 8736) , ld(50) 
common cl ,xl , lay.ntl ,ri ,ro,curc,kk,i i ,jj,nc 
format ( 1 x , 2 i 2 ) 

format (2f 10 .4,3i5,3f10.4,2i2) 

10 format (2f10. 4,2f12.10) 

11 format (3i4,6f5. 3) 

1 2 format (3i 3) 

18 format ( 3 i 3 , f 8.3) 



LUJ 


1 . 1 . ) 



6 


19 format (3f8. 2 ) 

111 format (1x,3f 6.2, 5x,3f 11 .2) ' ' . 

113 format (4x , 1 h X , 5x , 1 hY ,5x , 1 hZ , 1 5 x , 2hBX , 9x , 2hBY ,9x, 2hBZ) 

116 forma t (6x , 6n i nches , 20x , Shgauss ) 

155 format (/, 18a 4) . 

156 format ( 18a4) 

159 format (24x, 1 8a4) . 

192 f ormat (" 1 ", 1 40 , "system summary output") » 

1 93 format ( //t20 , "coi 1 descr i pt i on : " , 1 6 ( /t20 , 1 8a4)// ) 

450 read(1 , 10,end=420)da,ams ,xkt,xmu 
c 

c For calculating field inside any coil put the ip equal to the 

c number of tnat coil, 

c 

read ( 4 , 4 ) i co i 1 , i p 
xmp=xmu/(4. * 3. 141 5926) 

read( 1 , 1 56 , end=420 ) (con f g( 1c) , i c=1 , 1 8) 
c 

c Input spatial description: 

c maximum number of x , y , an d z i ncrements ,de 1 ta x, delta y.delte z 

c corner points of the plot, 

c 

read ( 1 , 1 1 ) i m , jm , km , dx , dy , dz , x ( 1 ) , y ( 1 ) , z ( 1 ) 
c 

c calculate and store the spatial coordinates for this calculation, 
do 200 0 i =2 , im 

2000 x( i )=x ( 1 )-( i -1 )*dx 
do 2001 j =2 , jm 

2001 y ( j ) =y ( 1 )-( j -1 )*dy 
do 50 j = 1 ,10 0 

50 sum( i , j ) =0 . 

read(1 , 12)ncoi 1 , inpop t.out 
wr> i te( 7,159) (confg( 1c) , i c=1 , 18 ) 
do 200 nc=1coi 1 ,ncoi 1 
write (6,299) nc 
c 

c Corresponding values for each 8 coil set should be read from one file, 
c 

i f (nc . gt .8) go to 28 
299 format (lx, "COIL NUMBER", i 1 0 ) 

read( 1 , 155) ( coi 1 ( i c ,nc) , ic=1 , 1 8) 
goto 29 

28 continue 
i f (nc . gt . 16) go to 27 
read ( 2 , 1 55 ) ( co i 1 ( i c , nc ) , ic=1 , 1 8) 
go to 29 

27 i f (nc, gt .24 ) goto26 

read (3 , 155) ( coi 1 ( i c,nc) , i c=1 ,18) 
go to 29 

26 read( 4 , 1 55 ) ( coi 1 ( i c , nc) , i c=1 , 1 8) 

read ( 4 , 5)c1 , x 1 , 1 ay, ntl, kk.ri.ro, cure, 1st epx , i stepr 
go to 31 0 

29 continue 
i f (nc . gt .8)go to 59 
read(1 , 18)mm ,nn,ni nt ,ccur 
go to6 1 

59 continue 

i f,(nc . gt . 16) go to 57 
read (2 ,18)mm,nn,nint,ccur 
go to 61 
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57 read( 3 , 1 8)mm , nn , n i nt . ccur 

61 continue 

do 00 j = 1 , n i nt 

do 80 i =1 ,4 

if(nc.gt.8) goto79 

read( 1 , 1 9) xx ( i , j ) , yy ( i , j ) , zz( i ,j) 

go to 80 

79 continue 

i f (nc . gt . 16) go to 77 

nead(2 , 1 9 ) xx ( i , j ) ; yy ( i , j ) , zz( i , j ) 

go to 80 

77 read( 3 , 1 9) x x ( i , j ) , yy ( i , j ) , zz ( i , j ) 

80 continue 
ntot=mm*nn*n i nt 

call coors( n i nt ,mm , nn , nt ot , xx , xl , x2 , 1 0) 
cal 1 coors( n i nt ,mm , nn , nt ot , yy , y 1 , y2 , 1 0) 
cal 1 coorsfn int,mm,nn,ntot,zz,z1 , z2, 10) 
c 

c If the Inside field is desired this part will be done defining nn2,mm2. 

c 

if(nc.ne.ip) go to 106 
im=0 

mm2=mm/2 

nn2=1 

do 99 n=nn2 , nn2 
do 99 ml = mm2 , mm, mm2 

do 99 i pp = 4 ,4 „ 

i f (ml -2)55, 55, 56 

55 m=1 

go to 58 

56 m=m1 t 

58 im= im+ 1 | 

ld( im) =4*nn* (m-1 )+4*n+ipp 

ldl=ld(im) 

x( im) = (xl ( Id 1 ) +x2 ( Idl ) )/2 

y(im)=(y1(ldl )+y2( Id! ) )/2 • , 

z( im)= (zl ( Idl )+z2( Idl ) )/2 
99 continue 

106 continue 
lm=ntot 

310 continue 

ns = 0 

do 1 99 i =1 , i m 

do 1 99 j = 1 , j m i , 

do 199 k=1,Km 

ns=ns+1 

write(6,298) ns 

298 format (lx, "NO. OF POINT", i 1 0 ) 

i f (nc. gt.24) go to 201 , 1 

bx=0 . 0 

by=0 . 0 

bz=0 . 0 

1 l*ld( i ) 

do 2 1 0 1 ■ 1 , 1 m 

i f (nc. ne. ip) go to 105 

1 f ( 1 .eq. ll)go to 198 

go to 105 

198 wr i te( 7, 104) i .xl ( 1 ) ,y1 ( 1 ) ,z1 ( 1 ) ,x2( 1 ) ,y2( 1 ) ,z2( 1 ) , 1 

104 format (2x,i2,2x,6(f7.3,2x),i7) 

go to 210 
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105 continue 

dk = 1 i 

a=(x1 ( 1 )-x( i ) )/39. 37 
b=(x2( 1 )-x( t ) )/39. 37 
c=(y1( 1 )-y( j ))/39.37 
d=(y2( V)-y( j ))/39.37 
e=(zl ( 1 )-z(k ) )/39.37 
f = (z2( 1 ) - 2 ( K ) ) / 39 . 37 
197 continue 

u=c*f-d*e ■! 

v=e*b-f*a 

w=a+d-b*c 

195 format ( 1 x , » I dd = " , i 3 , 2x ,6 ( f 7.3 , 2x , “ i =" , 1 3) , " 1 ■ • , 1 3,2x , •' i «« , i 3) 
r 1 ■ (a* a+c*c+e*e ) ** .5 

r2= ( b* b+d+d+ f * f ) * * . 5 ^ . ' f 

rs=r 1 + r2 
rm=r1+r2 
rdraa*b+c*d+e*f 
r xr* ( u **2+v* *2+w+ + 2 ) ** . 5 
194 continue 

h« ( rm+ rdr )/r m 
1 f (h-0.01 )2 , 1 , 1 

1 g«rs/( rm*(rm+rdr) ) 
goto 3 

2 g= ( ( rs ) * ( rm- rdr ) ) / ( rm*rx r *rxn) 

3 curmsxmp*ccur*g*1 0000*dk 
bx1=curm*u 

by 1 =cu nn* v „ . 

bz 1 =curm*w 

b::“‘oxrbx1 

by=by+by1 
bz=bz+bz1 
2 1 0 con t i n ue 

if(nc. ie.24) goto 202 

201 call c i rce( x ( i ) , y ( j ) , z ( k ) ,bx,by ,bz,br , i stepx , i stepr ) 

202 continue 

sum( 1 , ns ) = sum( 1 , ns ) +bx 

sum( 2 , ns ) =sum( 2 , ns )+by . 

sum(3 , ns ) =sum( 3 ,ns )+bz 
i f ( nc . 1 e .24 ) go to 203 
sum( 6 , ns ) =sum( 6 , ns )+bx 
sum( 7 , ns ) =sum( 7 , ns )+by 
sum(8,ns)=sum(8,ns )+bz 

203 cont i nue 

199 continue 

200 continue 

wr i te ( 7,192) 

write(7,159) (conf g( ioc) , ioc=1 , 18) 

write(7,193) ( ( co i 1 ( i c , n i ) , i c= 1 , 1 8 ) , n i = i co 1 1 , nco i 1 ) 
wr1te(7,113) 
wr 1 te( 7, 1 16) 
ns = 0 

do 41 0 i =1 , i m 
do 41 0 j = 1 , j m 
do 41 0 k=1 , km 
ns=ns+1 

write (7, 111) x(i),y(j),z(k) f (sum(ii,ns),ii=1,3) 

410 continue 
420 continue 
stop 
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end 

c 

c This subroutine calculates the field of a solenoids due to circulan 

c current elements. 

c i stepx , istepr are the number of wires in x and r direction which 
c approximated as one wire, 

c 

subrou tine circe(xp,yp,zp,bx,by,bz,br,i stepx , i stepr ) 
douole precision xk1,xk2,k,e 
dimension r( 100) ,xw(100) 

common cl , xl , lay, nt 1 ,ri ,ro,curc,kk,i i ,jj,nc 

p i =3*a tan( sqr t (3. ) ) 

cc=. 0000002 

a=nt 1 - 1 

do 10 i » 1 , n t 1 

b= i -1 

10 r( i )=( ri+b*( ro-ri )/a)/39.37 

aa=lay 
aa = x 1 / aa 
do 20 i =1 , 1 ay 
x i = I 

x i =x i - . 5 

20 xw ( i ) =aa*xi 

br = 0 
bx = 0 

do 30 1 * 1 , 1 ay , i stepx 
do 30 j = 1 ,nt 1 , istepr 

rho=( sq'rt (zp**2+yp**2) J/39.37 *. 

c~ { xw ( i J+cl ) *kk 
xx=(xp-c)/39 .37 
a = r( j ) 

xk2=4. *a*rho/( (a+rho) **2+xx**2 ) 

k=pi/2 

e=k 

i f ( xk2 . eq. 0 ) go to 25 
i f ( xk2 . eq. 1 ) go to 30 
xkl =1 -xk2 

g0 = 13. 86294361 1 2d- 1 
g1=0.96663442590d-1 
g2=3. 590092383 d-2 

g3 = 3. 7425637 13d-2 "* 

g4 = 14.511 962 1 2d~3 
g5= . 5d 0 

g6= 1 24 . 98593 597d-3 
g7=.6880248576d-1 
g8=332 . 83553 4 6d-4 
g9=441 .78701 2d-5 
h0 = 443 2 .5141 463d-4 
h1=.626060122d-1 
h2 = 4 . 7 573835 46d-2 
h3=1 .7 365064 51 d-2 
h4 = 24 . 998368 3 1 d-2 
h5=.9200180037d-1 
h6=40. 69697526d-3 
h7=. 5264496 3 9d~2 

k=g0+g 1 *xk1 +g2*xk1 **2+g3*xk1 * * 3+g4*xk1 **4+ ( g5+g6 *xk1 +g7*xk1 **2 
1 +g8*xk 1**3+g9*xk1**4)*(1og(1./xk1)) 
e=1+h0*xk1+h1*xk1 * *2 + h2* xkl **3+h3*xk1 **4+( h4*xk1 +h5*xk 1 **2+h6*Xk1 
1 **3+h7 *xk1 ** 4) * ( log( 1 ./xkl ) ) 

br=cc* curc*xx*1 000 0/( rho* ( (a+r ho) **2+xx**2 ) ** . 5) * (-k+{ a* *2+ 
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1rho**2+xx**2)*e/((a-rho) **2+xx **2) )+br 
25 continue 

bx=bx+cc*curc*10000/(( (a+nho) * *2+xx**2)»* . 5)*(k+(a**2-rho**2~ 

1 xx**2 ) *e/ ( (a-rho) * *2+xx* *2) ) 

30 continue 

xs=istepx 

rs=istepr " 

br=or* xs*rs 

bx=ox*xs+rs 

by=0 . -‘ ; 

bz = 0 

i f ( rho . eq .0 ) go to <10 
theta=atan(yp/zp) 
bz=bp*cos(theta) 
by=Dp* s i n( tneta) 

40 continue 

retunn 
knd 
c 

c This subrout ine calculates the endpoints of wire filaments 

c in a coil. The coil is appnoximated by a series of quadri 1 ateral 

c interfaces, 

c variable definition 

c ww(i,) vertex coordinates for each interface, 

c wl(ind) first endpoint of filament, 

c w2(ind) second endpoint of filament, 

c 1 im second dimension of ww array. 

subroutine coors( n i nt.,m, n , ntot , ww, wl , w2 , li m) 
dimension ww (4, 1 im) ,w1 (ntot ),w2(n tot) 
i f (m. g t . 1 . and . n .gt . 1 ) go to 100 
i nd=0 

i f (m+n ,ne.2) go to 110 
wr i te( 7, 200 ) 

200 format ("the filament coordinates are already specified") 
return 

110 if(m.ne.l) go to 120 
do 40 j = 1 , n 
do 40 k = 1 ,ni nt 
ind=ind+1 

wl ( i nd ) = ( ww ( 1 ,k)*(n-j)+(j-1)*ww(4,k))/float(n-1) 

40 continue 
qo to 130 
120 do 50 i =1 ,m 

do 50 k = 1 , ni nt 
i nd= i nd+1 

wl ( ind ) = (ww( 1 , k )* (m- i ) + ( i -1 ) *ww(2,k) )/f loa t (m-1 ) 

50 continue 

go to 130 
100 ind=0 

do 10 i=1 ,m 
do 10 j=1 ,n 
do 10 k = 1 ,ni nt 
i nd= 1 nd+1 

w1(ind)=(ww(1,k)*(n-j)*(m-i)+ww(4,k)*(j-1)*(m-i)+ 

1 ww(2,k)*(i-1)*(n-j ) +ww (3,k)*(i-1)*(j-1)) 

2 /f loat( (m-1 )*(n-1 ) ) 

10 continue 

c 

c w2 is generated from wl so that the endpoints for each coil 




0 






c are together, 
c 


130 k=m*n 

nn = n i n t-1 
do 20 1 = 1 , k 
11 = 11-1 )*nint 
w 2 (l*nint)=w 1 (ll+ 1 ) 

do 20. i=1 ,nn 
w2( 1+1 1 ) =w1 ( i+ 1 +l 1 ) 
20 continue 
return 
end 
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Side and Lift; Force Coils Magnetic Field 
system summary output 



<3 





coi 1 

descr i pt i on : 





lift 

coi 1 

x(+) y(-) 

z(~) 




1 i ft 

CO i 1 

x(+) y(+) 

Z( + ) 




lift 

coi 1 

x(+) y(-) 

z( + ) 




lift 

co i 1 

X ( + ) y ( +) 

z(-) 




1 i ft 

coi 1 

x ( 

-) y(+) 

z(-) 




1 i ft 

co i 1 

x ( 

-) y(-) 

z(-) 




1 i ft 

coi 1 

x ( 

-) y( +) 

Z( + ) 




lift 

coi 1 

x ( 

-) y(-) 

z( + ) 




lift 

co i 1 

2 

x(+) y( 

-) Z(-) 




1 i f t 

coi 1 

2 

x(+) y(+) z( + ) 




1 i ft 

coi 1 

'2 

x(+) y( 

-) z(+) 




lift 

coi 1 

2 

x(+) y(+) z(-) 

**• 



1 i ft 

coil 

2 

x(-) y(+) z(-) 




1 i ft 

coi 1 

2 

x(-) y( 

-) z(-) 




1 i ft 

coi 1 

2 

x ( — ) y(+) z(+) 




lift 

coi 1 

2 

x(-) y( 

-) z(+) 


X 

Y , 

Z 



BX 

BY 

3Z 

i ncnes 




gauss 



5.00 

0.00 

0.00 



0.00 

0.00 

3210.41 

4.00 

0.00 

0.00 



0.00 

0.00 

2498.65 

3.00 

0.00 

0.00 


- 

0.00 

0.00 

1827.49 

2.00 

0.00 

0.00 


- 

0.00 

0.00 

1194.03 

1 .00 

0.00 

0.00 



0.00 

-0.00 

589.29 

0.00 

0.00 

0.00 


- 

0.00 

-0.00 

0.00 

-1.00 

0.00 

0.00 


- 

0.00 

-0.00 

-589.29 


b c 
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Drag Coils Magnetic Field 
system summary output 


coil description: 
drag coil 1 
drag coi 1 2 


X 

Y 

z 

BX 

BY 

i nches 


gauss 


5.00 

0.00 

0.00 

-2986.28 

0.00 

4.00 

0.00 

0.00 

-2441 .84 

0.00 

3.00 

0.00 

0.00 

-1862.83 

0.00 

2.00 

0.00 

0.00 

-1257.09 

0.00 

1.00 

0.00 

0.00 

-633.15 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

1 .00 

0.00 

0.00 

633.15 

0.00 


3Z 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

o.go 









Magnetizing Coils Magnetic Field 
system summary output 


coil description: 
Magnetizing coil 1 
Magnetizing coi 1 2 


X 

Y 

z 

BX 

BY 

i nches 


gauss 


5.00 

0.00 

0.00 

-7258.84 

0.00 

4.00 

0.00 

0.00 

-7326.72 

0.00 

3.00 

0.00 

0.00 

-73 79 i 74 

0.00 

2.00 

0.0 0 

0.00 

-7417.71 

0.00 

i .00 

0.00 

0.00 

-7440.53 

0.00 

0.00 

0.00 

0.0 0 

-7448.14 

0.00 

1.00 

0.00 

0.00 

-7440.53 

0.00 


bz 

0 . 00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
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